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ABSTRACT 

The very large telescope (VLT) interferometer (VLTI) in its current operating state is equipped with high-order 
adaptive optics (MACAO) working in the visible spectrum. A low-order near-infrared wavefront sensor (IRIS) 
is available to measure non-common path tilt aberrations downstream the high-order deformable mirror. For 
the next generation of VLTI instrumentation, in particular for the designated GRAVITY instrument, we have 
examined various designs of a four channel high-order near-infrared wavefront sensor. Particular objectives of 
our study were the specification of the near-infrared detector in combination with a standard wavefront sensing 
system. In this paper we present the preliminary design of a Shack-Hartmann wavefront sensor operating in 
the near-infrared wavelength range, which is capable of measuring the wavefronts of four telescopes simultane- 
ously. We further present results of our design study, which aimed at providing a first instrumental concept for 
GRAVITY. 

Keywords: VLT interferometer, near-infrared wavefront sensing, Shack-Hartmann sensor, HAWAII detector 

1. INTRODUCTION 

GRAVITY^ is an interferometric imager with 10 /xas astrometric capability, which will coherently combine the 
light from all four ESO VLT 8m unit telescopes. GRAVITY will carry out the ultimate empirical test whether or 
not the Galactic Centre harbours a 4 million solar mass black hole, without resorting to theoretical assumptions. 
If the current interpretation of the near-infrared (NIR) flares from SgrA* is correct, GRAVITY has the potential 
of directly determining the space-time metric around this black hole. The instrument will also be able to 
unambiguously detect and measure the mass of black holes in massive star clusters throughout the Milky Way 
and in many active galactic nuclei to z^O.l. It will make unique measurements on gas jets in young stellar objects 
and active galactic nuclei. It will explore binary stars, exo-planet systems and young stellar disks. Because of its 
superb sensitivity GRAVITY will excel in milli-arcsecond phase-referenced imaging of faint celestial sources of any 
kind. Because of its outstanding astrometric capabilities, it will detect motions throughout the local Universe 
and perhaps beyond. Because of its spectroscopic and polarimetric capabilities it is capable of detecting gas 
motions and magnetic fleld structures on sub-milliarcsecond scales. 

GRAVITY will measure fringe visibility magnitudes and fringe positions with high precision in order to achieve 
its goals. The visibility magnitude signal-to-noise ratio depends, among others, on the optical quality of the 
wavefronts brought to interference. The seeing-limited optical quality can be enhanced with the help of adaptive 
optics. Adaptive optics requires sufliciently bright reference stars to work properly. The VLTI MACAO^ systems 
provide adaptive optics using the visible light of reference stars. For bright, yet highly embedded (reddened) 
targets, as found, e.g., in the galactic center region, wavefront sensing at near-infrared wavelengths is required as 
long as laser guide stars are not available. Thus GRAVITY will employ near-infrared wavefront sensors, which 
in combination with the existing MACAO deformable mirrors form a complete near-infrared sensitive adaptive 
optics system for the VLTI^ and its suite of instruments. 
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The selected baseline concept for the infrared wavefront sensor is a Shack Hartmann system. It will be located 
in the VLTI laboratory, thus also sensing tunnel seeing (see section [t]) in the VLTI optical train. The wavefront 
correction will be applied to the MACAO deformable mirrors located at the Coude focal station of each unit 
telescope. The main requirements for this near-infrared adaptive optics system are the delivery of a 36% Strehl 
point spread function in K-band while guiding on a K=7 star. We discuss the pros and cons of using a single 
detector for all four telescopes rather than one detector per telescope. The requirements for the near-infrared 
detector are discussed and compared with available devices. Since pupil stability is an important issue for the 
VLTI, we present a scheme to measure possible pupil drifts in real-time with the wavefront sensor. 
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Figure 1. Schematic overview of GRAVITY. 



2. ADAPTIVE OPTICS REQUIREMENTS FOR GRAVITY 

The requirements on the adaptive optics (AO) system resulting from the GRAVITY science cases can be sum- 
marized as fohows: 

• provide wavefront correction 

• provide NIR wavefront sensing 

• ahow for off- axis wavefront sensing 

• wavefront sensor wavelength range: H- and K-band 



Image quality in terms of Strehl ratio (SR) requirements as a function of reference star brightness in K-band 
assuming median seeing conditions on Paranal are: 

• irk = 7: SR=0.25 at 45 degrees zenith distance and 7 arcsec distance to the reference star 

• mK= 7: SR= 0.4 at 45 degrees zenith distance on- axis 

• irk = 10: SR=0.1 at zenith and on-axis 

A possible galactic center scenario for observations with GRAVITY is shown in fig. [2] 
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Figure 2. The GRAVITY galactic center observing field with guide stars. 



3. SYSTEM BASELINE 

3.1 Wavefront sensor options 

In our design study for GRAVITY we have selected as baseline system a Shack-Hartmann type wavefront sensors 
(SHWFS) because of their robustness and uncomplicated assembly. In addition, and unlike, e.g., curvature or 
pyramid sensors, they provide a direct image of the reference source, which provides a useful diagnostics both 
during the commissioning as well as the science operation of the NIR AO system. Pyramid WFS are in general 
optimized to focus the light from a single star on the tip of the pyramid. Thus they are probably less optimal for 
wavefront sensing in crowded regions, like, e.g., the centers of starburst clusters, or in the case of observations 
of (close) binary stars. 



3.2 Detector options 

The selection of the basehne for the NIR wavefront sensor (WFS) detector for GRAVITY was mostly based on 
availability, and of course technical requirements. The baseline design follows the line of injecting all 4 telescope 
beams into 1 single detector rather than into 4 separate ones. Table [T] shows the requirements of the IR WFS 
Detector on which the selection process was based. 



Parameter 


Desired Characteristics 


Minimum Requirements 


Operating Range 


1 - 2.5 m 


1.45 - 2.45 m 


Frame size 


4 X 100 X 100 


4 X 60 X 6 X 6 


Frame rate 


up to 500 fps 


100 fps 


Pixel clock 


up to 2.5 MHz 


110 kHz 


Read noise 


< 10 e- 




Linearity 


< 1% after flat-fielding 




Uniformity 


< 10% 





Table 1. Infrared wavefront sensor detector requirements. 



After an intense market survey, and contact with different manufacturers, the following options were inves- 
tigated in detail: the HAWAH 1 and the new HAWAH 1-RG and 2-RG devices from Teledyne. The HAWAH-1 
focal plane array was discarded at an early stage of the Phase A study because the desired frame rates (up to 
500 fps) would not have been achieved (the maximum frame rates obtained would be 25 Hz for the largest frame, 
and 115 Hz for the optimized frame). 



Detector Type 


Frame size (pix) 


Pixel clock (kHz) 


Frame rate (Hz) 


HAWAH-1 


4x 100x100 


250 


25 




4x 60x6x6 


250 


115 


HAWAH- IRG / 2RG 


4x 100x100 


2500 


500 




4x 100x100 


1000 


200 




4x 100x100 


500 


100 




4x 60x6x6 


1000 


926 




4x 60x6x6 


550 


509 




4x 60x6x6 


250 


230 




4x 60x6x6 


110 


102 



Table 2. Comparison HAWAU-l vs. HAWAH- IRG/ 2 RG. 



On the other hand, the new -RG devices have better quantum efficiencies^ and allow clocking frequencies up 
to 5MHz,^ what fulfils the frame rate requirements. Since their video output signals are organized into stripe-like 
channels that can be read out in parallel, and the clocking direction can also be set, an area selected between 
2 channels can be read as twice as fast. In this way, each of the 4 beams can be spread over 2 channels and 
the readout time can thus be reduced by a factor two. In order to minimize line skipping, the beams should be 
positioned either close to the top or bottom edge of the array. 

From the operating point of view of the GRAVITY WFS, the performance of the HAWAII- IRG and the 
HAWAII-2RG arrays is the same, and therefore both could be used. There are at least 8 Channels needed 
to accommodate the 4 beams and both detectors fulfil this requirement. Since the HAWAII-2RG is already 
supported by ESO, the baseline design is based on it. 

3.3 Deformable mirrors 

To complete a fully working adaptive optics system, our baseline design includes the usage of existing VLTI 
infrastructure, in particular the MACAO deformable mirrors. MACAO uses bimorph deformable mirrors with 
60 actuators. 
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Figure 3. HAWAII- IRG channel distribution with best beam location possibilities. 

4. PHASE A ADAPTIVE OPTICS PERFORMANCE SIMULATIONS 

4.1 Estimation of the optimal bandpass 

The choice of the bandpass has to take into account: the background spectrum, the object spectrum, the read- 
out noise of the detector (RON) ,the sampling frequency Fs. The bandpass is optimized by choosing (Ai,A2) 

that make the SNR maximum: f^^ 0{X)dX/yJ J^^ B{X)dX + Fs.RON'^, where 0(A) is the object spectrum per 
second, B{X) the background spectrum per second, Fs the sampling frequency and RON the detector read-out 
noise. 

In the Galactic Centre case, modelling IRS7 (see fig. [2| by blackbody with T=840 K, we have studied the 
optimal bandpass versus sampling frequency (100 to 200 Hz) and RON (from 5 to 20 e~) for different VLTI 
temperatures. We have assumed a VLTI transmission of 0.1 (ie a VLTI emission of 0.9). As a conclusion, (i) the 
J band is not absolutely required, (ii) a good compromise could be to use the band 1.44-2.45 /im, (iii) a blocking 
filter in the 1.80-1.96 /im region, where the atmosphere becomes opaque , is not required. 

4.2 Dimensioning the wavefront sensor 

We have considered IRS 7 as the AO reference source and the optimal bandpass defined above. Building on 
experience with NAOS, we have chosen a pixel size equal to the width of the diffraction pattern of a subaperture, 
i.e. (X/dsubap)' This is also equal to the number of subapertures times the telescope diffraction limit. We have aslo 
assumed the following parameters: Tvlti=0.1^ Twfs=^'8^^ Tfiiter=^'^^ QE=0.6, tyLT/=12°C (corresponding 
to a background=26500 e~/m^/s/arcsec^), wind speed=12 m/s corresponding to r=3ms in the visible. We have 
then computed the error budget (fitting, aliasing, temporal, noise, anisoplanatism errors) for different reference 
source magnitudes, different seeing conditions, different RON and F^, and different number of Shack Hartmann 
subpupils. Results are given in the following figures, expressed in rd^. The Strehl ratio is computed on the PSF 
image: it has been calculated by a Monte-Carlo code, simulating phase screens, compensated by the influence 
functions of the MACAO mirror, introducing noise (propagated though the system control matrix taking into 
account the Hartmann geometry), aliasing properties, servo-lag and anisoplanatism. 

In these simulations, the sampling frequency and loop bandwidth have been optimized in order to minimize 
the phase error. In particular, at low flux, the bandwidth has been reduced in order to smooth the noise: this 
results in a higher temporal error, although the sum of (noise + temporal) is minimized. 

The SR estimates differ slightly from exp(— a^) at low values because the expression SR = exp(— a^) is just an 
approximation and above all because the various error sources (fitting, aliasing, temporal, noise, anisoplanatism) 
are correlated between each other, and thus do not simply add up together. 



mK=6.55 


seeing 


0.93" 


RON = 10e-, 400 Hz 


RON = 20 e-, 200 Hz 


11x11 


10x10 


9x9 


8x8 


11x11 


10x10 


9x9 


8x8 


Fitting 


0.57 


0.57 


0.57 


0.57 


0.57 


0.57 


0.57 


0.57 


aliasing 


0.06 


0.08 


0.12 


0.22 


0.06 


0.08 


0.12 


0.22 


temporal 


0.07 


0.07 


0.07 


0.07 


0.09 


0.09 


0.09 


0.09 


Noise 


0.09 


0.08 


0.06 


0.06 


0.13 


0.11 


0.09 


0.09 


Aniso 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


exp(-a^) 


42.12 


41.75 


41.02 


37.30 


39.60 


39.61 


39.21 


35.38 


SR 


42.12 


41.75 


41.02 


37.30 


39.60 


39.61 


39.21 


35.38 



Table 3. Error budget for IRS 7 in median seeing conditions. 



mK=6.55 


seeing 


1.50" 


RON = 10 e-, 400 Hz 


RON = 20 e-, 200 Hz 


11x11 


10x10 


9x9 


8x8 


11x11 


10x10 


9x9 


8x8 


Fitting 


1.25 


1.25 


1.25 


1.25 


1.25 


1.25 


1.25 


1.25 


aliasing 


0.14 


0.18 


0.26 


0.48 


0.14 


0.18 


0.26 


0.48 


temporal 


0.18 


0.18 


0.18 


0.18 


0.22 


0.22 


0.22 


0.22 


Noise 


0.14 


0.11 


0.11 


0.12 


0.23 


0.19 


0.19 


0.22 


Aniso 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


exp(-o^) 


15.52 


15.21 


14.05 


11.11 


13.55 


13.44 


12.46 


9.71 


SR 


15.52 


15.21 


14.05 


11.11 


13.55 


13.44 


12.46 


9.71 



Table 4. Error budget for IRS 7 in bad seeing conditions. 



mK=8.5 


seeinc 


0.93" 


RON = 10 e-, 100 Hz 


RON = 20 e-, 50 Hz 


11x11 


10x10 


9x9 


8x8 


11x11 


10x10 


9x9 


8x8 


Fitting 


0.57 


0.57 


0.57 


0.57 


0.57 


0.57 


0.57 


0.57 


aliasing 


0.06 


0.08 


0.12 


0.22 


0.06 


0.08 


0.12 


0.22 


temporal 


0.6 


0.58 


0.56 


0.55 


1.01 


1.09 


0.9 


0.98 


Noise 


0.28 


0.24 


0.20 


0.17 


0.25 


0.26 


0.23 


0.22 


Aniso 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


exp(-o^) 


20.52 


21.35 


21.95 


20.57 


14.04 


12.57 


15.16 


12.73 


SR 


25 


26 


26 


25 


19 


19 


20 


17 



Table 5. Error budget for mK=8.5 in median seeing conditions. 



mK=8.5 


seeing 


1.50" 


RON = 10 e-, 100 Hz 


RON = 20 e-, 50 Hz 


11x11 


10x10 


9x9 


8x8 


11x11 


10x10 


9x9 


8x8 


Fitting 


1.25 


1.25 


1.25 


1.25 


1.25 


1.25 


1.25 


1.25 


aliasing 


0.14 


0.18 


0.26 


0.48 


0.14 


0.18 


0.26 


0.48 


temporal 


1.3 


1.3 


1.25 


1.2 


2.2 


1.92 


1.99 


1.9 


Noise 


0.34 


0.30 


0.27 


0.25 


0.34 


0.43 


0.34 


0.27 


Aniso 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


0.16 


exp(-a^) 


4.13 


4.11 


4.13 


3.53 


1.68 


1.94 


1.84 


1.72 


SR 


6.4 


6.4 


6.4 


5.5 


4.0 


4.0 


4.5 


3.3 



Table 6. Error budget for mK=8.5 in bad seeing conditions. 



mK=10.5 


seeing 0.93" 


RON = 10 e-, 25 Hz 


11x11 


10x10 


9x9 


8x8 


Fitting 


0.57 


0.57 


0.57 


0.57 


aliasing 


0.06 


0.08 


0.12 


0.22 


temporal 


5.46 


3.74 


2.56 


1.9 


Noise 


0.32 


0.42 


0.41 


0.54 


Aniso 


0.07 


0.07 


0.07 


0.07 


exp(-a^) 


0.15 


0.76 


2.41 


3.68 


SR 


3.6 


6.0 


7.5 


7.7 



Table 7. Error budget for mK=10.5 in median seeing conditions. 



From these simulations, we concluded that the top-level requirements are met for the mK=7 cases and closely 
met for the mK=10 case. This holds for the above mentioned configuration of 10 e- RON, sampling frequency and 
number of sub-apertures depending on brightness of the guide star. For the GC case the top level requirements 
are always fulfilled, even with a RON of 20 e-. 

5. OPTICAL DESIGN 

The optical design of the WFS follows the general concept of providing one single system for all four beams, 
using UTl-4 or ATl-4. These beams when leaving the switch yard are separated by 240mm, each pupil diameter 
is 18mm. 

The first component of the WFS optics (see Fig. [4| is the pupil-de-rotator, one individual component for 
all four telescopes. This device is realized by a carefully aligned K-mirror assembly, the individual mirrors are 
protected silver coated. 

The following fiat mirrors M4-1 to M4-4 together with M5-1 to M5-4 are feeding the four beams into the 
common optics. At the same time, M4-1 to M4-4 serve as field selection unit, the intermediate pupil image 
position has to be adjusted to these mirrors. 

The following Achromat Al is producing a common focused image centered to a cold field stop. This field 
stop reduces stray-light and ghost images down to a negligible value. A small achromat A2 re-images the pupil 
onto the lenslet array, which is mounted directly in front of the detector array. The current optical design of 
these two achromatic lenses is producing a pupil diameter of 2.4mm, which has to be adjusted to the finally 
chosen pixel size, number of sub-pupils and detector array channels. 




Figure 4. Optical layout of the WFS: the right sub-image shows details near the lenslet array. 



6. INSTRUMENT CONTROL ELECTRONICS 



6.1 Requirements 

The instrument control electronics (ICE) realizes all control and monitoring tasks for the GRAVITY instrument. 
It controls several electro-mechanisms, sensors, cryostats and electronics cabinets. GRAVITY uses a large number 
of electro- mechanical devices, which increases the complexity of the control electronics. In order to keep the 
electronics design simple and allow for easy maintenance, standard actuators, sensors and control devices from 
same model are used. Whenever possible, ESOs standard components are used. 

6.2 Architecture 

The ICE is divided into several subsystems, reflecting the functional subunits of GRAVITY. The control electron- 
ics architecture is based on the computer system architecture of ESO. The instrument workstation is connected 
to the control electronics via LAN sub-net. Inside the cabinet the CPU of a Local Control Unit (LCU) com- 
municates directly with other devices via VME-bus. Non VME devices communicate with the CPU board via 
RS232 interface. 

Three LCUs are present in the system. One LCU controls the AO system and the Switch-yard, the other 
two LCUs control the Fibre Coupler including the Field Selector, the Beam Combiner, the Spectrometers, and 
the Metrology system. The figure below shows the control system architecture for the AO system. 



Instrument 
Workstation 



GRAVITY 
ICE 



Real time 
Computer 
{SPART^ 



Detector 
Readout 
Electronics 
(NGC) 



LCU 



Interface 
Boards 



Interlock 
system 



Control 

and 
Monitor 
Devices 



AO- ICE 



WFS 
Detector 



Switchyard 
and WFS 

Sub-sytems 



GRAVITY 



Figure 5. Control electronics architecture of GRAVITY. 



7. VLTI INTERNAL SEEING AND ITS IMPACT ON THE LOCATION OF THE 

WAVEFRONT SENSOR 

In order to characterise the wavefront aberrations caused by turbulence in the VLTI tunnels and light ducts, 
a measurements campaign was started in the last week of August 2007. The setup consists of an commercial 
Shack-Hartmann (SH) sensor (Wavescope, manufactured by Adaptive Optics Associates) with control PC and 
an additional frame grabber PC capable of recording the spot pattern of the SH sensor at a frame rate of 30 Hz. 
The control PC of the SH sensor is used for alignment of the Wavescope unit, checking of illumination level 
and integration times, and to decompose the recorded wavefront measurements into Zernike modes. While this 
system allows data acquisition only at a maximum speed of 15 Hz, the external frame grabber solution directly 
records the video signal of the Wavescopes internal camera at a full frame rate of 30 Hz. The recorded spot 



patterns can be fed back to the analysis software of the Wavescope unit to allow the same data reduction at 
doubled temporal sampling. For the example measurement discussed below, the HeNe laser beacon at the Coude 
focal station of the VLT was used as light source. The used delay line number 2 was set-up to compensate an 
optical path length of 100 m. In order to match the beam size (80 mm) to the entrance pupil size of the SH 
sensor (10 mm), a single converging lens was used. The SH sensor was equipped with a microlens array with 
480 jam pitch size, resulting in useful data in ^250-300 sub-apertures. 

7.1 Spatial wavefront characteristics 

In case of Kolmogorov turbulence, the variance of the wavefront (t| is related to the aperture size D, the Fried 
parameter ro, and the sum of the Zernike coefficient variances by: 



1.0299 



5/3 



X] < (tti - ai)2 > 



This assumes an infinite outer scale Lq, and it has to be decided from case to case if this assumption is valid 
for the given conditions. If the outer scale is finite but still well larger than the Fried parameter ro, the above 
formula can still be applied by using only the higher order Zernike modes, in the simplest case starting with the 
focus term. In this case the factor 1.0299 reduces to 0.134, reflecting that the major part of wavefront variance 
is contained in the tip/tilt modes. By comparing the two values for rO derived in this way, one can already 
test if Kolmogorov statistics with infinite outer scale or the van-Karman model with finite Lq allows the best 
representation of the measured turbulence. 

The following figure shows the temporal evolution of the wavefront rms with and without tip-tilt aberrations 
during a 5000 s measurement with 30 Hz sampling rate. Static aberrations have been removed. 



RMS wavefront error vs time 



RMS wavefront error without tip/tilt vs time 




Figure 6. Wavefront rms with static aberrations excluded (left) and without tip-tilt (right). 

The root mean square (rms) value of the wavefront is 126 nm, equivalent to 1.25 rad at the laser's wavelength 
or 0.36 rad in iT-band. The resulting Strehl ratio - assuming an otherwise perfect optical system - can be 
approximated by SR = exp"^'^^. 

The numbers for this example are SR=0.21 and SR=0.88 in the red and in i^-band, respectively. At the laser 
wavelength of 632.8 nm, the corresponding Fried parameter ro for this series is 6.3 cm. Taking the beam compres- 
sion ratio of ^100 into account, this translates to an effective on-sky ro of 6.3 m. At 2.2/im, the corresponding 
value would be 28.1 m (using the relation ro ~ A^/^), well larger than the telescope diameter. 

If a perfect tip/tilt-correction can be accomplished, the remaining wavefront rms can be evaluated by considering 
only the higher Zernike modes, starting with the focus term (see right plot of Fig. |6|. In this case, the mean 
wavefront rms reduces to 43.9 nm, giving an effective on-sky value for ro of 29 m in the red and 130 m in K-band. 



The fact that the wavefront rms after removal of tip and tilt was reduced by a factor of 2.87 shows that the 
turbulence is quite similar to Kolmogorov turbulence where a factor of 2.8 is expected (variance reduces from 
1.0299 to 0.134). 

7.2 Temporal wavefront characteristics 

Figure [7] shows the temporal power spectrum of the total wavefront rms evolution. 

^ ^ Temporal PSD of the rms wavefront error 
Temporal PSD of the rms wavefront error ^ q2 | , , , , , , , , , ^ , , , , 




0.001 0.010 0.100 1.000 10.000 100.000 5 10 15 

Frequency [Hz] Frequency [Hz] 

Figure 7. Temporal power spectrum of the wavefront RMS. The left plot uses logarithmic scaling on both axes, the right 
one only for the power. 

8. SUMMARY 

Adding near-infrared wavefront sensing capabilities to the VLTI will increase the number of possible science 
targets for VLTI instruments using adaptive optics. While compensating high-order atmospheric turbulence up 
to the Coude focal station of the UTs (location of the deformable mirror) there are currently two options available 
to compensate residual low-order aberrations up to the VLTI laboratory with the IRIS wavefront sensor. One 
uses the XY table of the MACAO wavefront sensor^ as an indirectly and rather slow tilt compensating device, 
the other option uses fast piezo driven steering mirrors inside the VLTI laboratory instrument feeding optics 
(currently only available for AMBER and PRIMA). For GRAVITY we investigated the power of high order 
aberrations introduced between the Coude deformable mirror (M8) and the interferometric lab. With typically 
around 50 nm rms it is questionable whether they need to be compensated. For the residual low-order aberrations, 
a compensation as close as possible to the science instrument seems mandatory in order to stabilize the beams 
within fibers and other position sensitive devices. 
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